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Abstract — Partial Unit Memory (PUM) codes are a special class 
of convolutional codes, which are often constructed by means 
of block codes. Decoding of PUM codes may take advantage 
of existing decoders for the block code. The Dettmar-Sorger 
algorithm is an efficient decoding algorithm for PUM codes, but 
allows only low code rates. The same restriction holds for several 
known PUM code constructions. In this paper, an arbitrary- 
rate construction, the analysis of its distance parameters and a 
generalized decoding algorithm for PUM codes of arbitrary rate 
are provided. The correctness of the algorithm is proven and it 
is shown that its complexity is cubic in the length. 

Index Terms — Convolutional codes, Partial Unit Memory 
Codes, Bounded Minimum Distance Decoding 

I. Introduction 

The algebraic description and the distance calculation of 
convolutional codes is often difficult. By means of block 
codes, special convolutional codes of memory m — 1 can 
be constructed, which enable the estimation of the distance 
parameters. Moreover, the existing efficient block decoders 
can be taken into account in order to decode the convolutional 
code. There are constructions of these so-called Partial Unit 
Memory (PUM) codes yj, ||2J based on Reed-Solomon (RS) 
||3)-||5), BCH ||6), Q and - in rank metric - GabiduHn |8|, 
Q codes. Decoding of these PUM codes uses the algebraic 
structure of the underlying RS, BCH or Gabidulin codes. 

In Dettmar and Sorger constructed low-rate PUM 

codes and decoded them up to half the extended row distance. 
Such a decoder is called Bounded Minimum Distance (BMD) 
decoder for convolutional codes. Winter pT| gave first ideas 
of an arbitrary rate construction. 

In this contribution, we construct PUM codes of arbi- 
trary rate, prove their distance properties and generalize the 
Dettmar-Sorger algorithm to PUM codes of arbitrary rate. We 
prove the correctness of the decoding algorithm and show that 
the complexity is cubic in the length. To our knowledge, no 
other construction and efficient decoding of PUM codes of 
arbitrary rate exist. Due to space limitations, we consider only 
PUM codes, but all results apply also to Unit Memory codes. 

This paper is organized as follows. In Section III] we 



give basic definitions. Section III provides the arbitrary rate 
construction and calculates its parameters. In Section |IV] we 
explain and prove the BMD decoding algorithm. Section |V] 
concludes this contribution. 

This work was supported by the German Research Council "Deutsche 
Forschungsgemeinschaft" (DFG) under Grant No. Bo 867/21-1. 



II. Definitions and Notations 

Let q be a power of a prime and let F denote the finite field 
of order q. We denote by F" = F^^" the set of all row vectors 
of length n over F and the elements of a vector a.j E F" by 



Let us define a zero-forced terminated convolutional code C 
for some integer L by the following Lk x (n{L+m)) generator 
matrix G over the finite field F 



G 



/ Go Gi . . . G„ 
Go Gi 



G. 



Go Gi 



(1) 



where G^, z = 0, . . . , m are k x n-matrices and m denotes 

dof 



the memory of C as in [12]. In the following, N = L + m. 

The error-correcting capability of convolutional codes is 
determined by extended (or active) distances. 

Let C^(j) denote the set of all codewords corresponding to 
paths in the minimal code trellis that diverge from the zero 
state at depth and return to the zero state for the first time 
at depth j. The extended row distance of order j is defined as 
the minimum Hamming weight of all codewords in C^{j): 



{wt(c)}. 



Similarly, let C^(j) denote the set of all codewords leaving the 
zero state at depth and ending in any state at depth j and let 
C'^^{j) denote the set of all codewords starting in any state at 
depth and ending in the zero state in depth j, both without 
zero states in between. The extended column distance and the 
extended reverse column distance are: 



d'v min {wt(c)}, d'-^ =^ min {wt(c)}. 



The free distance is the minimum (Hamming) weight of 
any non-zero codeword of C and can be determined by 
dfioc — minjidj}. The extended row distance can be lower 
bounded by a linear function with slope a: 



a = lim (-^|. 

i-i-oo L j J 



PUM codes are convolutional codes of memory to = 1. 
Therefore, the semi-infinite generator matrix consists of two 
k X n sub-matrices Go and Gi. Both matrices have full rank 



2 



Go = 



(2) 



if we construct an {n,k) UM code. For an {n,k \ ki) PUM 
code, rank(Go) = k and rank(Gi) — ki < k hold, such that: 

Goo \ t~< ^ 1 

Go J ' ^'-[ 

where Goo ™d Gio are ki x n matrices and Gqi is a (fc — 
ki ) X n-matrix. The encoding rule for a code block of length 
n is given by Cj — ij ■ Gq + • Gi, for ij, ij^i £ F*^. 

The free distance of UM codes is upper bounded by di^cc < 
2n — k + 1 and of PUM codes by d^cc < n — k + ki + 1. For 
both the slope is upper bounded hy a < n — k Q, (Ts). 

As notation, let the generator matrices 



Goi 
Gio 



Goi and G^ = 




define the block codes Co, Ci, Coi and Ca with the minimum 
Hamming distances do, di, doi and da and the BMD block 
decoders BMD(Co), BMD(Ci), BMD(Coi) and BMD(Ca), 
which correct errors up to half their minimum distance. 

III. Constructing PUM Codes of Arbitrary Rate 
A. Construction 

Since each code block of length n of the PUM code can 
be seen as a codeword of the block code Ca, a great da is 
important for the distance parameters of the convolutional code 
as well as for the decoding capability. One approach is to 
define by Gq a Maximum Distance Separable (MDS) code and 
da = n—k—ki + 1. This is basically the construction from |j6), 
ITpl which designs low-rate PUM codes since the (fc + fci) x n 
matrix Gq can define an MDS code only if fc + fci < n. 
Otherwise (as observed by |11 1), there are linear dependencies 
between the rows of Gq, what we have to consider when 
constructing PUM codes of arbitrary rate. In the following, 
we provide a construction of arbitrary fci < fc and calculate 
its distance parameters. 

Let fc + fci — < n, for some ip < ki, and let the (fc + fci — 
ip) X n matrix 



Goi 



with the sub-sizes 



A : {ki — (p) X n 
^ : ip X n 
Goi : (fc — fci) X n 
B : {ki — If) X n 



(3) 



define an MDS (e.g. RS) code. We define the sub-matrices of 
the semi-infinite generator matrix of the PUM code as follows 
in order to enable arbitrary code rates. 

Definition 1 (PUM Code of Arbitrary Rate) Let fci < fc < 

n and let Gtot be defined as in ([3]). Then, we define the PUM 
code by the following submatrices (|2|.- 



Go 



Goo 
Goi 



Gi = 







(4) 



Since Gtot defines an MDS code, Co, Ci and Cio (compare 
Section [n] for the notations) are also MDS codes. We restrict 
tp <ki since otherwise all rows in Gi are rows of Gq. Note 
that any rate k/n in combination with any fci is feasible with 
this restriction since fc + l<fc + fci — (/j<n and hence, we 
have only the trivial restriction k < n. 



B. Calculation of Distances 

We calculate the extended row distance of the construction 
from Definition [T] by cutting the semi-infinite generator matrix 
into parts. Each code block of length n can be seen as a 
codeword of Cq with minimum distance 



d{Ga) = d{Gtot) 



fc — fci + (/3 + 1. 



However, due to the linear dependencies between the sub- 
generator matrices, a non-zero information block can result in 
a zero code block. The following lemma bounds the maximum 
number of such consecutive zero code blocks. 

Lemma 1 (Consecutive Zero Code Blocks) The maximum 
number i of zero code blocks Cj, Cj+i, .... Cj+£_i, which have 
no edge in common with the zero state, is 



ki- (p 



Proof: If ip = 0, there is no zero code block obtained 
from a non-zero information block and £ — 0. 

For < ip < ki, let 



-i,0,...,0|«fei 




((y.. j), zo,...,zfc,_^_i , |o^^^) 

(£-l)(fci-v) fci-v ip-{l-l)(ki-ip) k-ki 

(0^^^^,io, . . • , V-£(fci-i^)-i I 0^^^^). 

k~ki 



ip-£{ki~ip) 



In the non-binary case, each second block ij, ij+2, ■ ■ ■ has to 
be multiplied by — L Then, 



= ■ Gi + ij+h • Go = 0, yh = 0, 



In each step, we shift the information vector to the right by 
ki — (p positions, where this shift size is determined by the 
size of A. Since $ has p rows, this right- shifting can be done 
\(p/{ki — (p)~\ times. We ceil the fraction since the last block 
ij+e-i can contain less than ki — ip information symbols. ■ 
Therefore, after £ zero code blocks there is at least one block 
of weight da and the slope can be lower bounded by: 



a > 



n — k — ki+p> + l 



(5) 



The extended distances can be estimated as follows. 



Theorem 1 (Extended Distances) The extended distances of 
order j for the PUM code of Definition [7] are: 

d'[>d[^ doi, >d'j ^ do + {j - 2) ■ a + di, j > 1, 
> d] = do + {j - 1) • a, j> 0, 

with doi — n — k + ki + 1, do = di — n^k + 1 and a as in (fSl) 
and dj, dj and d^ denote the designed extended distances. 



3 



Proof: For the calculation of the extended row distance, 
we start in the zero state, hence, the previous information 
is Iq = 0. We obtain dT^ for an information block ii = 
(0, . . . , 0, i]^^ , . . . , j^lli), then Ci e Cqi . The extended row 



distance of order j follows from (|5]l and a last information 



block ii 



(0,. 



i'-''^ 



U) 



The second-last block 



is arbitrary and thus Cj — ij ■ Gq + ij-i ■ Gi is in Ci. 

The calculation of the extended column distance starts in 
the zero state, hence, Iq ~ 0, but we end in any state, thus, 
di > do- For higher orders, each other block is in Cq. 

The reverse extended column distances considers all code 
blocks starting in any state, hence there is no restriction on 
io,ii and Ci £ Ca- In order to end in the zero state, ij — 
(0, . . . , 0, i]^^ , . . . , i^^li) and as for the extended row distance 
Cj e Ci. ■ 

The free distance is then the minimum, i.e., 

rfficc ^ niiii {^i} = min{n — k + ki + 1,2 ■ (n — k + 1)}. 

i— 1,2,... 

Note that if dficc = ri — fc + fci + 1, then the free distance is 
optimal since the upper bound is achieved Q. 

IV. BMD Decoding Algorithm 

A. BMD Condition and Idea 

Let the received sequence r = c+e = (rp, ri, . . . , r^^i) be 
given, where r/j = C/j + e/j, h = 0, . . . , N — 1 is in F", c = 
(cg, Ci, . . . , Cjv-i) is a codeword of the (terminated) PUM 
code as in Definition [T] and is an error block of Hamming 
weight wt(eft ). A BMD decoder for convolutional codes is 
defined as follows. 

Definition 2 (BMD Decoder for Convolutional Codes ||10|) 

A BMD decoder for convolutional codes guarantees to find 
the Maximum Likelihood (ML) path as long as 



E 



wt(e/i) < 



(6) 



holds for all j = 0, 



, iV — 1 and i — 1, 



Algorithm [T] shows the basic principle of our generalization 
of the Dettmar-Sorger algorithm to arbitrary rate. 

Algorithm 1: Arbitrary-Rate Decoder for PUM codes 
Input: Received sequence r of length N ■ n 

1 Decode block Tq with BMD(Co), 

decode blocks Tj for j = 1, . . . , W - 2 with BMD(Ca), 
decode block rjv_i with BMD(Ci), 
calculate ij if £ + I consecutive blocks were decoded 
successfully and assign metric as in ([T] 

2 From all found blocks ij, decode steps forwards with 
BMD (Co) and steps backwards with BMD(Ci) 

3 From all found blocks ij, decode next block with 
BMD(Coi) and assign metric as in 

4 Search the complete path of smallest weight with the 
Viterbi algorithm 

Output: Information sequence i of length {N — 1) ■ k 

The main idea of the algorithm is to take advantage of the 
efficient BMD block decoders for Ca, Cq, Ci and Cqi. With 



the results of the block decoders, we build a reduced trellis 
and finally use the Viterbi algorithm to find the ML path. 
Since this trellis has only very few edges, the overall decoding 
complexity is only cubic in the length. Figure [T] illustrates the 
decoding principle for £ — 1. 

r ~1 

Step 1 

y I I I * * I I I 



step 2 




Fig. 1. Example of the decoding algorithm for £ = 1. where the three first 
steps of Algorithm ^ for the received sequence r are illustrated. 

Since each code block of the PUM code of length n 
is a codeword of the block code C^, the first step of the 
algorithm is decoding with BMD(Cq.). Due to the termination, 
the first and the last block can be decoded with BMD (Co), 
respectively BMD(Ci). The decoding result of BMD(Ca) is 
Cj. Assume it is correct, then Cj = Cj 



j — («q"' ^\ . . . ^i'^^^^) is a part of the previous 
information block. Now, we want to reconstruct the informa- 
tion ij — («o''\ ■ ■ ■ :*fc'-i) ™d ^^"^ '■^i^' need 
consecutive decoded code blocks since the linear dependencies 
"spread" to the next I blocks as shown in Example [T] 



where i 



1-1 ] 



1 



Example 1 (Reconstructing the Information) Let 

Lp = 2/3ki, where £ ~ 2 and $ has twice as much 
rows as A. Assume, we have decoded Co, Ci and C2 
and we want to reconstruct ii. Decompose io,ii,i2 into: 
ij = (iW I il^l I if I if) for j = 0, 1, 2, where the first three 



sub-blocks have length ki 



(f and the last k — ki. Then, 



ct - (ifl 



where $ 



.[2] , .[1] I .[3] , .[2] I .[4] I .[3]^ 



dcf^ 



*2 

Goi 

) and $1, $2 have ki — if rows. Since we 
know Ci and Gtot defines an MDS code, we can reconstruct 
the vector ii. This directly gives us if and if. This can 
be done in the same way for Cq and we also directly obtain 
(among others) if. To obtain if, we substract if from the 
known sum if +ig^'. For C2, this reconstruction provides if 



and we have the whole ii. This principle also gives us i. 



[ki 



[1] I ;[21 I .[3] 



ip ). This is why £+1 consecutive decoded blocks 



are necessary to reconstruct an information block. Note that it 
does not matter if the other decoded blocks precede or succeed 
the wanted information, this principle works the same way. 



4 



After this decoding and reconstruction, we build an edge in 
a reduced trellis for each block with the metric: 



wt(rj — Cj) if Step 1 finds Cj and ij, 
L(da + 1)/2J else. 



(7) 



Remark 1 The error of minimum weight causing a sequence 
of non-reconstructed information blocks in Step 1 is as follows: 



(0,...,0,x I 0,...,0,x I ... I 



£+1 blocks 



i+1 Mocks 



0,.-.^,0,x I 0,...,0), 

l+l blocks I blocks 



where the x marks blocks with at least dajl errors. Also the 
information of the error-free blocks cannot be reconstructed, 
since we need l+l consecutive decoded blocks. The last I 
error-free blocks are the reason why we substract £ in the 
definitions of and L^pK This corresponds to i additional 
decoding steps in forward direction. The (minimum) average 
weight in a sequence of non-reconstructed information blocks 
(without the last £ blocks) is therefore da/{2{£+ 1)). 

Assume, in Step 1, we decoded c, and reconstructed and 
a part of the previous information then we calculate: 



Gn 



' 



■ ' 'fel-1 ) 



Gn 



•(j-l)|-(j-2) 



(8) 



Hence, as a second step, we decode Ip'' blocks forward with 
BMD(Co) respectively blocks backward in BMD(Ci). 
These codes have higher minimum distances than da and close 
(most of) the gaps between two sequences of correctly decoded 
blocks in Cq. The values £p^ and £^^^ are defined by: 



£^i^^ = mill 



da 



h=l 



mm 

i=l,2,.. 



da 



rrij+h ^ d^ 
- 2 

z 



£ + 1 



ij-h 



h=l 



£ + 1 



> 



(9) 



(10) 



Lemma |3] in Section [IV-B| proves that after Step 2, the size of 
the gap between two correctly reconstructed blocks is at most 
one block. 



For Step 3, assume we know i^^^^j^ 



from Step 1 and ij_2 from Step 1 or 2, then similar to ([8]l 



which shows that we can use BMD(Coi) to close the remaining 
gap at j — 1. After Step 3, assign as metric to each edge 



L(doi + 1)/2J 



if BMD(Co), BMD(Ci) or BMD(Co 



else, 



(11) 

where again Cj denotes the result of a successful decoding. 
Note that there can be more than one edge in the reduced 
trellis at depth j. 



Finally, we use the Viterbi algorithm to search the ML path 
in this reduced trellis. As in |10|, we use nij as edge metric 
and the sum over different edges as path metric. 

Section IV-B proves that if (|6]l is fulfilled, after Steps 1-3, 
all gaps are closed and Algorithm [T] finds the ML path. 
It is a generalization of the Dettmar-Sorger algorithm to 
arbitrary rates, which results in linear dependencies between 
the submatrices of the PUM code (see Definition [T}. This 
requires several non-trivial modifications of the algorithm. 
Namely these are: the reconstruction of the information 
requires £ + 1 consecutive code blocks (see Example [T]l, 
the path extensions (|9]l, ( [TOj i have to be prolonged and the 
assigned metric has to be adapted appropriately (|7]i, ( fTTj ) since 
the smallest error causing a non-reconstructable sequence is 
generalized as in Remark [T] 

B. Proof of Correctness 

In this subsection, we prove that Algorithm [T] finds the 
ML path if (|6]) is fulfilled. For this purpose. Lemma [2] shows 
that the size of the gaps after Step 1 is not too big and in 
Lemma [3] we prove that after Step 2, the gap size is at most 
one block. Finally, Theorem |2] shows that we can close this 
gap and that the ML path is in the reduced trellis. Then, the 
Viterbi algorithm will find it. The complexity of the decoding 
algorithm is stated in Theorem |4] 

Lemma 2 The length of any gap between two correct 



reconstructions 

min(4\4+^)) 



'« Step 1, ij, 
if ^ holds, with 



is less than 



mm 

1=1,2,.. 



mm 

1=1,2,.. 



E 

i 

E 

h=l 



da — TTljj^h 

£ + 1 

da — nij- 
£ + 1 



> 



> 



Proof: Step 1 fails if there occur at least da/ 2 errors in 
every {£ + l)-th block, followed by £ correct ones (compare 
Remark Assume there is a gap of at least Lp'^ blocks after 
Step 1. Then, 



^wi{eh) > ^ 2( 

h=l h=l ^ 



1) - ^ £+1 - 2 

' h=l 



contradicting (|6]l. We prove this similarly for L^g without 
substracting I in the limit of the sum, since we directly start 
left of the £ correct blocks on the right. Therefore, the gap 
size is less than min(L^\ i^''). ■ 

Lemma 3 Let ij and ij+i be reconstructed in Step 1. Let 



Step 2 decode £p^ blocks in forward and £g^^' blocks in 
backward direction (see (|9]l, jlO) ). Then, except for at most 
one block, the ML path is in the reduced trellis if ^ holds. 

Proof: First, we prove that the ML path is in the reduced 
trellis if (|6| holds and in each block less than min{(io/2, di/2} 
errors occurred. In this case, BMD(Co) and BMD(Ci) will 
always yield the correct decision. The ML path is in the 
reduced ti-ellis if t,^' + ' > i - 1, since the gap is then 

< i — 1 and at least da/ 2 



closed. Assume that £ 



5 



errors occur in every {£+ l)-th block in the gap, since Step 1 
was not successful (compare Remark [T]l. Then, 

i-l 



> 



h=l 



dfU) 
2 

do 



dnU + i) 



')da 



+ (z - 3) 



da 



2(^+1) 



2{£ + l) 



which is a contradiction to (|6]). 

Second, we prove that at most one error block e^j, j < h < 
j + i has weight at least (io/2 or di/2. To fail in Step 1, there 
are at least da/ 2 errors in every {£ + l)-th block. If two error 
blocks have weight at least do/ 2 = di/2, then 



^wt (e,+,)>2.^ 

h=l 



da 

2 



> 



z 



in contradiction to (|6]l. Thus, the ML path is in the reduced 
trellis except for a gap of one block. ■ 

Theorem 2 (|6]l holds, the ML path is in the reduced trellis. 

Proof: Lemma [3] guarantees that after Step 2, the gap 
length is at most one block. This gap can be closed in Step 3 
with Coi, which is always able to find the correct solution 
since doi > di — dftce- ■ 



C. Decoding of a Single Block 

Similar to pO) , we give a weaker BMD condition to 
guarantee ML decoding of a single block. This condition 
shows how fast the algorithm returns to the ML path after 
a sequence where (|6]) is not fulfilled. A BMD decoder for 
convolutional codes guarantees the correct decoding of a block 
Yj of a received sequence r = c + e if the error e satisfies 

k+i-l -^r 

2-1. (12) 



E 

h=k 



wt (eh) < y , yi,k with k<j<j 



To guarantee ( [T2] l for a certain block if ^ is not fulfilled for 
the whole sequence, we introduce an erasure node in each step 
j as in |7|, representing all nodes which are not in the reduced 
trellis. Let ej, ej^i denote erasure nodes at time j, j — 1 and let 
Sj, Sj-i be nodes found by BMD decoding in Steps 1 and 2. 
Let tB denote the minimum number of errors of any edge 
starting from Sj_i and Sj in forward, respectively backward 
direction, ta denotes the minimum number errors of any edge 
between nodes at time j — 1 and j. We set the metric of the 
connections with the erasure nodes as follows. 



Connect 


Metric 






= misj-i) + 


max([{do + l)/2J, do-tp) 






m{sj) 


= miej-i) + 


max(L(di+l)/2J, dt-tg) 






m{ej) 


= ^{€^-1) + 








{da ~ ta) 

, L(rf<. + 1)/2_ 


if 3 an edge between Sj-i,Sj 
, else. 



Theorem 3 If ( [T2| i holds for rj, the Viterbi algorithm for the 
reduced trellis with erasure nodes finds the correct block Cj. 



Proof: The metric of the erasure nodes is always at least 
/2. All nodes of a state are connected with the erasure nodes 
of the previous and the next state. As soon as ( [T2| ) is fulfilled, 
the metric of a correct edge is better than all other edges and 
the ML path will be chosen. ■ 

D. Complexity Analysis 

The complexity is determined by the complexity of the 
BMD block decoders, which are all in the order 0{n^), if 
the construction is based on RS codes of length n. 

Similar as Dettmar and Sorger we can give the 

following bound on the complexity. Due to space restrictions, 
the proof is omitted here. 

Theorem 4 Let C be a PUM code as in Definition [7] where 
Gtot is the generator matrix of an RS code. Then, the decoding 
complexity of Algorithm^of one block is upper bounded by 

CpuM <o{ie + i)dan^) ^o{ie + l)n^) . 

V. Conclusion 

We presented a construction of PUM codes of arbitrary rate 
and provided and proved an efficient decoding algorithm. The 
algorithm corrects all error patterns up to half the designed 
extended row distance, where the complexity is cubic in the 
length of a block. For £ — 0, the Dettmar-Sorger algorithm 
1 10 1 is a special case of Algorithm [T] 
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